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£NJ ■ Abstract. The use of orthogonal polynomials for integral models on the lattice is applied 

p. ' to the 3nj-symbols that appear in the coupling of several angular momenta. These symbols 

, are connected to the Ponzano-Regge method to solve the Einstein equations on a discrete 



Riemannian manifold. 



> ■ 1. Classical orthogonal polynomials of one discrete variable 

These polynomials satisfy a difference equation of hypergeometric type such that the 
C\l ■ difference derivatives of the some polynomial satisfy a similar equation. The discrete variable 
can be consider of two types. 

a) On homogeneous lattice: x = 0,1,2, ■■ ■ 
45 The corresponding polynomials satisfy a difference equation 

Oh' 

45 

B 

> 



a(x)AVp n (x) + r(x)Ap n (x) + X n p n {x) = 

withA/(x) = f(x+l) — f(x), Vf{x) = f(x) — f(x — l), a(x) and r(x) are functions 
of second and first order respectively; 
an orthogonality relation 
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with p(x)a weight function and d n a normalization constant. To these polynomials 
correspond the Meixner, Kravchuk, Charlier and Hahn polynomials [1] 

b) On non homogeneous lattice: x = x(s), s = 0,1, 2, • • • the polynomials satisfy a 
difference equation 

A Vy n (x) 1 f Ay n (x) Vy n (x)\ 

Ax s-h ] Vx(s) 2 I Ax{s) Vx (s I 



Yl y™( S ) Vm{ s ) Pn(x) Ax ( S - ]-) =d\ 5 nm . 



2. 

an orthogonalit relation 
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v„(s) v^(s) n„(x) Ax ' ,s — 2 

The corresponding polynomials are classified according to the latice function: 
For x(s) = s(s + 1), we have the Racach and dual Hahn polynomials 
For x(s) = q s or we have the g-Kravchuck, g-Meixner, g-Charlier and g-Hahn 

polynomials. 

For x{s) = g3+ 2 g or 2 — ^ - , we have the g-Racach and g-dual Hahn polynomials [2]. 
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2. Generalized Clebsch-Gordon coeficients and generalized 3nj-Wigner symbols 

If two angular momentum operators are coupled to give a total angular momentum J = J1 + J2 
the new basis can be expressed in terms of the old ones 

\j1j2jm) = (jihm 1 m 2 I 31323m) \jij2m 1 m 2 ) . 

mi+m2=m 

The symmetry properties of the Clebsch-Gordon coefficients in this expansion are more patent 
if one sbstitutes them by the Wigner symbols 

(3i32jm I j 1 j2m 1 m 2 ) = (-l) jl - j2+j - m yj2jTi ^ J m ) 

Similarly if we couple three angular momentum operator we obtain a new basis in terms of 
the old ones: 

\j1j2j3j12jm) = J2 (ji32j3m 1 m 2 m 3 | jij 2 j 3 ji2jm) \jij 2 j3m 1 m 2 m 3 ) 
for the coupling ( J\ + J 2 ) + J3 = J, 

\j1j2j3j23jm) = (ji3233m 1 m 2 m 2 , \ 3132333233m) \3i3233m 1 m 2 m 3 ) 

for the coupling Ji + ( J 2 + J 3 ) = J. 

Both bases are related by some matrix U (J12, J23) that can be written in terms of 
generalized 6j-Wigner symbol 



U(j 12 ,j23) = {-l) h+h+j * +j J(2h2 + 1) (2J23 + 1) 



32 312 
33 3 323 



In similar fashion can be written the generalized Clebsch-Gordon coefficients and generalized 
3nj -Wigner symbols [3]. The algebraic properties of these symbols can be represented by 
geometrical graphs [3]. 

3. 3nj-symbols as orthogonal polynomials of several discrete variable 

The 6j- symbols are proportional to the Racah polynomials through the following relation [2] 



(-l) Jl+J2+J2 V(2ji 2 + l)(2j 2 3 + l){(| I g ) } = ^^\x, a, b) 
with 

X(s) = S(S + 1) , S= 323 

a = 33 ~ 32 , b = 3 + j 3 + 1 , n = j 12 - ji + j 2 - j 
a = ji- 32 -js+j , P = ji - 32 + 33 ~ j 

Using the assymptotic limit of the Racah polynomials and the connections between the 
Jacobi polynomials and the Wigner little functions one can prove the following approximation 
of the 6j-sumbols when ji ~ j 2 ~ j 3 ~ j jn 

f 3i 32 312 1 _ I L ) d n2 ( ^ (l) 

1 33 3 323 J - VJ1+J2 + W33 + 3 + 1 n - n ' XA ' 3 
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with 

a (2jg3 + I)' ~ (jl + 32 + I)' ~ (jg + j + I)' 
C °S v = — : — — : — 

2(jl+j2 + l)(j3+J + l) 

The 3nj-symbols of the first and second kind can be written in terms of 6j-symbols, and 
therefore in terms of product of Racah polynomials, giving rise to orthogonal polynomials of 
several discrete variables. To illustrate this take, f.L, the 12j-symbol of the second kind as a 
combination of 6j-symbols. 



ji h 
k k 
h k 2 
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Here R n = {ji + k + h). Substituting each 6j-symbol for the corresponding Racah 

i=i 

polynomial we obtain: 



f ji h is U } i 4 JpOA . a . 

\k k h k \= £ tt—t II (k) = Pn (k k k k) 

Ul h k 3 kj * ZX + L i=l d r H 

which is a polynomial of four discrete variables. 
For the assymptotic limit we find 

f ji h is j4 ) 4 i 

Ii k k k } ^V(2x + l)TT- ; -d%_ k .:_ k u 

1 , , , , I — ' V ' A. I U. I 1 J* K nJt + l K i + 1( mod 4) v ' 

\ ki k 2 k 3 k A ) * %=i l 

These formulas can be easily generalized to any 3nj-symbols of first and second kind. 



E(2* + l)(-l 



,«„ + !.( J Ji h x 
k 2 ji k 



J2 ^2 

h J3 



/ 2 



J3 ^3 ^ 

^4 J4 ^3 



ji h x 
h ji k 



4. Application to spin networks and to Ponzano-Regge integral action 

Penrose has proposed a model for the space and time in which the underlying structure is 
given by a set of interactions between elementary units that satisfy the coupling of angular 
momentum operators, called spin networks [4]. One particular case of these networks can be 
described by the graphs of 3nj-symbols. From different point of view Regge has proposed a 
method to calculate Einstein action by the approximation of curved riemannian manifold by a 
polyedron built up of triangles. Later Ponzano and Regge applied the properties of 6j-symbols 
to calculate the sum action over this triangulation [5] 

Let M be a riemannian manifold that is approximated by a polyedron with boundury D 
and it is decomposed into p tetrahedra T k represented by 6j-symbols. 

The polyedron give rise to triangular faces /, represented by 3 j-symbols, and to q internal 
edges xu as well as to external ones U with respect to the boundary D. 

Ponzano and Regge define the sum 

s=En r *(- i rn( 2a *+ i ) ^ 

Xi k = l 1 = 1 

When li — > oo, % — > 0, Uli — > finite we recovered the continuous manifold. In order to 
compute the 6j-symbols in te classical limit, we uses the assymptotic formula [2] 
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<Ln> (*) « ("I) 



/ 2 / 2j + m - m' + 1 N 
7r (j — m) \ 2j — m + ml + 1 



cos 



(;+*)«-( 



m — m 



Vsin (9 



at m ~ m' ~ 1, j » 1. Substituting this expression in (1) with m — j± — j 2 , m' = 
J4 — is 5 J — J6 and taking the edges of the tetrahedra j x + \, ■ ■ ■ , 3% + \, very large except 



j 6 , we have 



Jl 32 33 



k 35 k) VT2nV IV 2 J V 2)2 

1 \ 7T / 1\7T / 1\7T 7T 



+ (J 2 +2J2-1 J4+ 2J2 + ^ 5 ' 27 2 ' 4 



6 



cos 



.1=1 



1\. 7T 



where 9i is the dihedral angle for the edge ji and V = | (ji + |J (j4 + |j (^6 + 1) sen ^> 
Note the formula (3) has been proved rigurously by Roberts [5]. Introducing formula (3) in 
formula (2), Ponzano and Regge proved that it leads in the continuous limit to the integral 
action of the general relativity. 
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